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Abstract 

We consider modeling for strong-strong beam-beam interactions beyond preceding lin- 
earized/perturbative methods such as soft gaussian approximation or FMM (HFMM) etc. 
In our approach discrete coherent modes, discovered before, and possible incoherent oscil- 
lations appear as a result of multiresolution/multiscale fast convergent decomposition in the 
bases of high-localized exact nonlinear modes represented by wavelets or wavelet packets 
functions. The constructed solutions represent the full multiscale spectrum in all internal hid- 
den scales from slow to fast oscillating eigenmodes. Underlying variational method provides 
algebraical control of the spectrum. 
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Abstract 

We consider modeling for strong-strong beam-beam in- 
teractions beyond preceding linearized/perturbative meth- 
ods such as soft gaussian approximation or FMM (HFMM) 
etc. In our approach discrete coherent modes, discov- 
ered before, and possible incoherent oscillations appear as 
a result of multiresolution/multiscale fast convergent de- 
composition in the bases of high-localized exact nonlinear 
modes represented by wavelets or wavelet packets func- 
tions. The constructed solutions represent the full mul- 
tiscale spectrum in all internal hidden scales from slow 
to fast oscillating eigenmodes. Underlying variational 
method provides algebraical control of the spectrum. 

1 INTRODUCTION 

We consider the first steps of analysis of beam-beam in- 
teractions in some collective model approach. It is well- 
known that neither direct PIC modeling nor soft-gaussian 
approximation provide reasonable resolution of computing 
time/noise problems and understanding of underlying com- 
plex nonlinear dynamics [1], [2]. Recent analysis, based as 
on numerical simulation as on modeling, demonstrates that 
presence of coherent modes inside the spectrum leads to 
oscillations and growth of beam transverse size and defor- 
mations of beam shape. This leads to strong limitations for 
operation of LHC. Additional problems appear as a result 
of continuum spectrum of incoherent oscillations in each 
beam. The strong-strong collisions of two beams also lead 
to variation of transverse size. According to [2] it is reason- 
able to find nonperturbative solutions at least in the impor- 
tant particular cases. Our approach based on wavelet anal- 
ysis technique is in some sense the direct generalization 
of Fast Multipole Method (FMM) and related approaches 
(HFMM). After set-up based on Vlasov-like model (ac- 
cording [2]) in part 2, we consider variational-wavelet ap- 
proach [3]-[17] in framework of powerful technique based 
on Fast Wavelet Transform (FWT) operator representa- 
tions [18] in section 3. As a result we have multiresolu- 
tion/multiscale fast convergent decomposition in the bases 
of high-localized exact nonlinear eigenmodes represented 
by wavelets or wavelet packets functions. The constructed 
solutions represent the full multiscale spectrum in all inter- 
nal hidden scales from slow to fast oscillating eigenmodes. 
Underlying variational method provides algebraical control 
of the spectrum. 
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2 VLASOV MODEL FOR BEAM-BEAM 
INTERACTIONS 

Vlasov-like equations describing evolution of the phase 
space distributions tj) 1 = ip 1 * (x,p x ,6) (j = 1,2) for each 
beam are [2] : 
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where 



p*(x,8)= / i/j*(x,p X) 0)dp x 



(2) 



and if>* is the density of the opposite beam, q x is unper- 
turbed fractional tune, £ x is horizontal beam-beam param- 
eter, N is a number of particles, x, p x are normalized 
variables. This model describes horizontal oscillations of 
flat beams with one bunch per beam, one interaction point, 
equal energy, population and optics for both beams. 

3 FWT BASED VARIATIONAL 
APPROACH 

One of the key points of wavelet approach demonstrates 
that for a large class of operators wavelets are good ap- 
proximation for true eigenvectors and the corresponding 
matrices are almost diagonal. FWT [18] gives the maxi- 
mum sparse form of operators under consideration (1). It 
is true also in case of our Vlasov-like system of equations 
(1). We have both differential and integral operators inside. 
So, let us denote our (integral/differential) operator from 
equations (1) as T (L 2 (R n ) -> L 2 (R n )) and its kernel as 
K. We have the following representation: 



<Tf,g>= 



K(x,y)f{y)g(x)dxdy (3) 



In case when / and g are wavelets tpj k — 2^ 2 ip{2^x — k) 
(3) provides the standard representation for operator T. Let 
us consider multiresolution representation ... C Vi C 



Vl C Vq C V-! C V-2 



The basis in each V< is 



tfij,k (x), where indices k, j represent translations and scal- 
ing respectively. Let Pj : L 2 (R n ) — > Vj (j € Z) be pro- 
jection operators on the subspace Vj corresponding to level 
j of resolution: (Pjf)(x) = X)fc < f>fj,k > fj,k( x )- Let 
Qj = Pj-i — Pj be the projection operator on the subspace 
Wj (Vj-i = Vj © Wj) then we have the following "micro- 
scopic or telescopic" representation of operator T which 



2 



takes into account contributions from each level of reso- 
lution from different scales starting with the coarsest and 
ending to the finest scales [18]: 

t = E (Qi T Qi + Qi TP i + P o T Qo)- w 

We remember that this is a result of presence of affine 
group inside this construction. The non-standard form of 
operator representation [18] is a representation of operator 
T as a chain of triples T = {Aj,Bj,Tj}j e z, acting on 
the subspaces Vj and Wf Aj : Wj — > Wj,Bj : Vj — > 
Wj,Tj : Wj — > Vj, where operators {Aj,Bj,Tj}j eZ 
are defined as Aj = Q 3 TQ V B } = QjTPj, Tj = 
PjTQj. The operator T admits a recursive definition via 

where Tj = PjTPj and Tj acts on Vj : Vj -> V$. It should 
be noted that operator A, describes interaction on the scale 
j independently from other scales, operators Bj , Tj de- 
scribe interaction between the scale j and all coarser scales, 
the operator Tj is an "averaged" version of Tj_\. We 
may compute such non-standard representations for differ- 
ent operators (including Calderon-Zygmund or pseudodif- 
ferential). As in case of differential operator as in other 
cases we need only to solve the system of linear algebraical 
equations. The action of integral operator in equations (1) 
we may consider as a Hilbert transform 

(Hp*)(x) = -p.v. ^ ^^-dx' (6) 

The representation of H on Vo is defined by the coefficients 

r e = J <p(x - £){H<p)(x)dx, I e Z. (7) 

which according to FWT technique define also all other 
coefficients of the nonstandard representation. So we have 
H = {Aj, Bj,Tj}j e z with the corresponding matrix ele- 
ments a,i-e, bi-t, Ci-t which can be computed from coeffi- 
cients re only: 

L-l 

a% = E 9kgk'r2i+k-k' 

k,k'=0 
L-l 

bi = ^2 9kh k 'T 2 i+k-k' (8) 

k,k'=0 
L-l 

Ci = E h k g k >r 2 i + k-k' 

k,k'=0 

The coefficients re (7) can be obtained from 

L/2 

rt =r 2 e + ^2 d 2k -i(r 2 e-2k+i + r 2e+2k -i) (9) 

k=l 



where d n are the so called autocorrelation coefficients 
of the corresponding quadratic mirror filter {hk} k Zo : 
d n = 2 Ylf=o~ n hih i+n , n = 1, . . . , L - 1, d 2k = 0, 
k = l,...,L/2-l,g k = {-l) k h L - k - 1 ,k = 0,...,L-l, 
which parametrizes the basic refinement equation <p(x) = 
y/2 J2k=a h k tp(2x — k). This equation really generates all 
wavelet zoo. It is useful to add to the system (9) the follow- 
ing asymptotic condition re = —l/n£ + 0{(r 2M ), which 
simplifies the solution procedure. Then finally we have 
the following action of operator Tj on sufficiently smooth 
function /: 

(Tjf)(x) = ]T [2-'' ]T nfj, k -A yj, k {x), (10) 

k£Z \ i / 
in the wavelet basis tpj_ k (x) = 2~ : >/ 2 tp(2~ : >x — k) where 

/,> : 2 ' j f(x)<p(2-ix-k + e)dx (11) 

are wavelet coefficients. So, we have simple linear para- 
metrization of matrix representation of our operator (6) in 
wavelet bases and of the action of this operator on arbitrary 
vector in proper functional space. The similar approach 
can be applied to other operators in (1). Then we may ap- 
ply our variational approach from [3]-[17]. Let L be an 
arbitrary (non) linear (differential/integral) operator corre- 
sponds to the system (1) with matrix dimension d, which 
acts on some set of functions ^> = ^(9, x,p x ) = 

(*\9,x,p x ), . . . ,^ d {e,x,p x )y e,x, Px e!lc R 3 , 

= L(Q,9,x,p x )^(9,x,p x ) = 0, where Q = 
Qd 1 ,d 2 ,d 3 (°i x ,Px,d/d9, d/dx,d/dp x ,J dxdp x ). Let us 
consider now the N mode approximation for solution as the 
following ansatz (in the same way we may consider differ- 
ent ansatzes) [17]: 

N 

y N (9,x, Px ) = a rsk A r ®B s ®C k (6,x,p x ) (12) 

r,s.k— 1 

We shall determine the coefficients of expansion from 
the following conditions (different related variational ap- 
proaches are considered in [3]-[16]): 

C m = J (L^ N )A k (9)Bt(x)C m (p x )d9dxdp x = 

So, we have exactly dN 3 algebraical equations for dN 3 
unknowns a rsk . The solution is parametrized by solutions 
of two set of reduced algebraical problems, one is linear 
or nonlinear (depends on the structure of operator L) and 
the rest are some linear problems related to computation 
of coefficients of algebraic equations. These coefficients 
can be found by some wavelet methods by using com- 
pactly supported wavelet basis functions for expansions 
(12). We may consider also different types of wavelets 
including general wavelet packets. The constructed solu- 
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tion has the following multiscale/multiresolution decom- 
position via nonlinear high-localized eigenmodes 

1>(e,x, Pai )= Y a ijk A i (6)B^x)C k (p x ), (13) 
A\6) = A% slow (6) + Y A i(u r &), u r ~ 2 r 

r>N 



Bi{x) = Br i l ° w {x) + Y Bj(kh), kf ~ 2 l 

1>M 

C s ( Px ) = Cr l ° w ( Px ) + Y C6 m(.k 2 mPx ), k 2 rn ~ 2" 



which corresponds to the full multiresolution expansion 
in all underlying time/space scales. Formula (U3b gives 
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Figure 1: Region of nonlinear resonances. 
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Figure 2: Eigenmodes decomposition. 

us expansion into the slow part fj}°M L and fast oscillating 
parts for arbitrary N, M, L. So, we may move from coarse 
scales of resolution to the finest one to obtain more detailed 
information about our dynamical process. The first terms 
in the RHS of formulae (13) correspond on the global level 
of function space decomposition to resolution space and 
the second ones to detail space. The using of wavelet basis 
with high-localized properties provides fast convergence of 
constructed decomposition (13). In contrast with different 
approaches, formulae (13) does not use perturbation tech- 
nique or linearization procedures and represents dynamics 



via generalized nonlinear localized eigenmodes expansion. 
Numerical calculations are based on compactly supported 
wavelets and related wavelet families. Figures 1,2 demon- 
strate resonances region and corresponding nonlinear co- 
herent eigenmodes decomposition according to representa- 
tion (13). 

4 REFERENCES 

[1] K.Yokoya, Phys.Rev. ST AB, 3, 124401, 2000, M.P. Zorza- 
no, F. Zimmerman, Phys.Rev. ST AB, 3, 044401, 2000 

[2] M.P. Zorzano, LHC Report 440, CERN, 2000 

[3] A.N. Fedorova and M.G. Zeitlin, Math, and Comp. in Simu- 
lation, 46, 527, 1998. 

[4] A.N. Fedorova and M.G. Zeitlin, New Applications of Non- 
linear and Chaotic Dynamics in Mechanics, 31, 101 Klu- 
wer, 1998. 

[5] A.N. Fedorova and M.G. Zeitlin, CP405, 87, American 
Institute of Physics, 1997. Los Alamos preprint, physics/- 
9710035. 

[6] A.N. Fedorova, M.G. Zeitlin and Z. Parsa, Proc. PAC97 2, 
1502, 1505, 1508, APS/IEEE, 1998. 



L7J 



|8] 



[9] 



[10] 



[12] 



[14] 



[15] 



A.N. Fedorova, M.G. Zeitlin and Z. Parsa, Proc. EPAC98, 
930, 933, Institute of Physics, 1998. 

A.N. Fedorova, M.G. Zeitlin and Z. Parsa, CP468, 48, 
American Institute of Physics, 1999. Los Alamos preprint, 
physics/990262. 

A.N. Fedorova, M.G. Zeitlin and Z. Parsa, CP468, 69, 
American Institute of Physics, 1999. Los Alamos preprint, 
physics/990263. 

A.N. Fedorova and M.G. Zeitlin, Proc. PAC99, 1614, 1617, 
1620, 2900, 2903, 2906, 2909, 2912, APS/IEEE, New York, 



1999. Los Alamos preprints: physics/9904039, 9904040 



9904041, 9904042, 9904043, 9904045, 9904046, 9904047. 

[11] A.N. Fedorova and M.G. Zeitlin, The Physics of High 
Brightness Beams, 235, World Scientific, 2000. Los Alamos 



preprint: physics/0003095 



A.N. Fedorova and M.G. Zeitlin, Proc. EPAC00, 415, 
1101, 1190, 1339, 232 5, Austrian Acad. Sci.,2000. 
Los Alamos preprints: ^hysics/0008045| , 



872, 



0008046, 



0008047, 0008048, 0008049, 0008050. 

[13] A.N. Fedorova, M.G. Zeitlin, Proc. 20 International Linac 
Con f., 300, 303, SLA C, Stanford, 2000. Los Alamos prepri- 
0008200. 



nts: physics/0008043 



A.N. Fedorova, M.G. Zeitlin, Los Alamos preprints: phy- 
sics/0101006, 0101007 and World Scientific, in press. 

A.N. Fedorova, M.G. Zeitlin, Proc. PAC2001, Chicago, 
790-1792, 1805-1807, 1808-1810, 1811-1813, 1814-1816, 
2982-2984, 3006-3008, IEEE, 2002 or arXiv preprints: 



3hysics/0 106022, 0106010, 0106009, 0106008, 0106007, 



0106006, 0106005. 

[16] A.N. Fedorova, M.G. Zeitlin, Proc. in Applied Mathematics 
and Mechanics, 1, Issue 1, pp. 399-400, 432-433, Wiley- 
VCH, 2002. 

[17] A.N. Fedorova, M.G. Zeitlin, this Proc. 

[18] G. Beylkin, R. Coifman, V. Rokhlin, CPAM,44, 141, 1991 



4 



